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Abstract 

The leading relativistic and recoil corrections to bound state ^r- factors of particles with arbitrary 

f~^ ■ spin are calculated. It is shown that these corrections are universal for any spin and depend only 

CN , on the free particle gyromagnetic ratios. To prove this universality we develop nonrelativistic 

W)! quantum electrodynamics (NRQED) for charged particles with an arbitrary spin. The coefficients 

in the NRQED Hamiltonian for higher spin particles are determined only by the requirements of 

Lorentz invariance and local charge conservation in the respective relativistic theory. For spin one 

charged particles the NRQED Hamiltonian follows from the renormalizable QED of the charged 

vector bosons. We show that universality of the leading relativistic and recoil corrections can be 

r~ I ! explained with the help of the Bargmann-Michael-Telegdi equation. 
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Gyromagnetic ratios of particles in hydrogenlike bound states have become in the last 
ten-fifteen years an active field of experimental and theoretical research. The gyromagnetic 
ratio of a bound electron is proportional to the ratio of the spin fiip and cyclotron frequencies 
of a hydrogenlike ion and to the electron-ion mass ratio. The experimental uncertainties of 
the ratio of the spin fiip and cyclotron frequencies of the hydrogenlike carbon ^'^C^+ and 
oxygen ^^O^"*" were reduced to 5 — 7 parts in 10^'', see [l|, l2| and review in |3|. The theoretical 
expression for the bound state (^-factor was also greatly improved recently (see, e.g. J4|,[5| and 
references in [3|), and the theoretical uncertainty was reduced to 1.5 — 5.5 parts in 10^^. As a 
result measurements of the bound electron (/-factor became the best source for precise values 
of the electron mass in atomic units [3|. This bright picture is marred by the discrepancy 
on the magnitude of the leading relativistic and recoil corrections to bound state f^-factors 
existing in the literature |6|-l8|. This discrepancy shifts the theoretical value of the bound 
state 5f-factors of the hydrogenlike carbon i2(^5+ ^^^ oxygen ^®0^+ by about 2 — 3 parts 
in 10^^. It will become even more phenomenologically relevant if proposed improvement 
[9| of the experimental accuracy by two orders of magnitude is achieved. Theoretically, 
discrepancy between different results for the leading relativistic and recoil corrections to 
bound state (/-factors is connected with different treatments of the spin dependence of these 
corrections. Below we derive an effective NRQED Hamiltonian for charged particles with 
arbitrary spins and calculate the leading relativistic and recoil corrections to the bound 
state (yf-factors in loosely bound two-particle systems. We show that these corrections are 
universal for all spins; they do not depend on the magnitude of spin. 

A loosely bound two-particle system is effectively nonrelativistic, with characteristic ve- 
locities of constituents of order Za. We are looking for the leading binding and recoil cor- 
rections of order (Za)^. NRQED is a natural framework for calculation of these corrections. 
We first consider leading nonrecoil corrections to the gyromagnetic ratio of order [ZaY- To 
calculate all such corrections we need the NRQED Lagrangian that includes all terms in 
nonrelativistic expansion up to and including f^. We should also include in the effective 
Lagrangian terms with the external Coulomb field Aq, since for such field (eAn) ~ (Za)^. 
The NRQED Lagrangian for spin one half case is well known (see, e.g., [lO|). The coef- 
ficients in the NRQED Lagrangian for charged particles with an arbitrary spin should be 
determined from comparison with the results of the respective relativistic theory. The prob- 
lem is that renormalizable QED for charged particles with high spin does not exist. The 
rules for calculation of all one-photon interactions of charged particles with arbitrary spin 



were constructed some time ago in [111 . |12| | . This construction uses only Lorentz invariance 



and local current conservation, and it should be valid for charged particles of arbitrary spin. 



The interaction vertex in the approach of [ll|, |12| is a direct generalization of the ordinary 
spin one half vertex 

where q = P2 — Pi, ^fiu is the generalization of ordinary spin one half auu, ^"^ is a covariant 
spin four-vector, r = (g ■ S)"^, and Fe(0,0) = 1, Fm(0,0) = g/2. The wave functions are 
spinors with dotted and undotted indices that are symmetrized among themselves (for more 



details see llHlSJ). The form of the vertex in Eq. ([T]) is uniquely fixed by the requirements 
of Lorentz invariance, C, P and charge conservation. Charged particles with higher spins 
automatically carry higher multipole moments that arise as coefficients in expansion of the 
form factors F^ and F^ over r. These intrinsic electric and magnetic multipole moments are 
treated phenomenologically, and we do not try to calculate them. The phenomenological 



approach to multipole moments is an advantage for our purposes because we would like to 
describe how g-iactois of not necessarily electromagnetic origin (for example the g-iactor of 
a spin one deuteron) change in a loosely bound electrodynamic system. 

In the spin one half case the NRQED Lagrangian is constructed from the gauge invariant 
operators Z) = V — ieA = i{p — eA), E, B, and S. For higher spin particles, besides 
the spin operator, we should also include higher irreducible intrinsic multipole moments as 
the building blocks of the NRQED Lagrangian. Technically these multipole moments are 
polynomials in the components of the spin operator that for higher spins do not reduce to 
numerical tensors and operators linear in spin. The most general NRQED Lagrangian has 



the form (compare with |10| ) 
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where Qij = SiSj + SjSi — {2/3)S'^Sij is proportional to the electric quadrupole moment 
operator {Qij = for spin one half), and is a 25* + 1-component spinor field for a particle 
with spin s. We included in the Lagrangian in Eq. ([2]) operators of dimensions not higher 
than four, except those (like the terms with derivatives of magnetic field) that are irrelevant 
for calculation of the leading recoil corrections. Let us mention that gauge invariant bilinears 
in E and B are of too high order to generate leading relativistic contributions of order {ZaY 
to bound state ^f-factors. 

The coefficients in Eq. ([2]) are usually determined from comparison of the one- and two- 
photon scattering amplitudes in NRQED and relativistic QED. Although some terms in 
Eq. 02]) are bilinear in A and E all such terms can be restored from one-photon terms 
due to gauge invariance. Then the one-photon relativistic vertex in Eq. ([1]) is sufficient 
for calculation of all the coefficients in Eq. (|2]). We calculated scattering amplitudes off an 
external electromagnetic field using the nonrelativistic Lagrangian in Eq. (|2]) and using the 
relativistic one-photon vertex in Eq. ([1]) at r = 0. In the relativistic calculation we used 
noncovariantly normalized particle spinors in the generalized standard representation, which 
is necessary for consistency with the respective nonrelativistic results. Diagrammatically this 
choice of spinors and representation corresponds to the Foldy-Wouthuysen transformation 
(for more details, see, e.g., 13[). After nonrelativistic expansion we compared results of the 



relativistic calculation with the nonrelativistic ones and obtained values of all constants in 
the Lagrangian in Eq. ([2]) 

q S^ 
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where S^ = 45*, A = 1/(25-1) for integer spin and S^ = 45" + 1, A = 1/(25) for half integer 
spin. Dependence on the magnitude of charged particle spin arose in the coefficients before 



the Darwin term and the induced electric quadrupole interaction. The (^-factor in Eq. ([!]), 
Eq. ([3]) is the total gyromagnetic ratio of a free nonrelativistic particle defined through the 
effective interaction Hamiltonian Hint = —ge/{2m)B ■ S. If the charged particle is subject 
only to electromagnetic interaction then g reduces to a sum of the QED perturbation series. 
For spin one half the coefficients in Eq. ([3j) coincide with the respective coefficients in [lO|, 
if the phenomeno logical gf-f actor is substituted in the expressions in [lOJ] instead of the per- 
turbative g = 2(1 + a/27T). As an independent test of the effective Lagrangian in Eq. ([2]) we 
considered the charged PF^-boson sector of the Glashow-Weinberg-Salam Electroweak The- 
ory amended by the anomalous magnetic moment term. We derived the effective NRQED 
Lagrangian for the W^ bosons. This Lagrangian coincides with the Lagrangian in Eq. ([2]) 
for spin one charged particles. 

The NRQED coefficients in Eq. ([3]) are calculated ignoring all relativistic loop diagrams 
and q^ and r dependence of the form factors in Eq. ([1]). Both the loop diagrams in relativistic 
QED and multipole expansion of the form factors would generate further corrections to the 
coefficients in Eq. (^. However, we are interested only in corrections to bound state g- 
factors of order (Za)"^ ~ v"^. The coefficients in the effective Lagrangian are calculated 
comparing scattering amplitudes in relativistic and nonrelativistic theories. Counting of 
powers of the coupling constants in the case of scattering amplitudes calculated at a generic 
kinematical point is trivial. In ordinary renormalizable spin one half QED (as well as in the 
renormalizable QED of spin one W^ vector bosons) all diagrams, besides those that give 
contributions only to the free particle ^f-factors, generate corrections to the coefficients that 
are additionally suppressed by powers of Za. We expect the same effect in any reasonable 
theory for higher spin particles. It is also obvious that accounting for g^ and r dependent 
terms in the form factors in Eq. ([T]) generates terms suppressed by additional powers of 
Za. We do not need to consider two-photon Compton effect diagrams, since all terms in 
Eq. ([2]) bilinear in fields can be restored from one-photon diagrams with the help of gauge 
invariance. Any gauge invariant terms connected with the two-photon diagrams are of too 
high order in Za to contribute to the leading relativistic corrections of order (Za)^. Hence, 
the Lagrangian in Eq. ([2]) with the coefficients from Eq. ([3]) is sufficient for calculation of 
the leading relativistic corrections to the bound ^f-factor in the nonrecoil case. 

Our goal is also to calculate recoil corrections of order [Za)"^ that are linear and quadratic 
in the mass ratio. To this end we need to construct an effective two-particle NRQED 
Hamiltonian for a loosely bound electrodynamic system of two particles. The interaction 
between two charged particles with accuracy up to [ZaY is described by the one photon 
exchange which generates Coulomb and Breit interactions. We calculated the one-photon 
potential for two particles with arbitrary spins and magnetic moments and obtained 
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where ?"i(2), Pi(2), 'S'i(2), ""^1(2), fi'i(2), and Q^^^- are the coordinate, momentum, spin, mass, 
gyromagnetic ratio, and induced quadrupole moment of the first (second) particle, and 
r = ri — r2 is the relative coordinate. 

This interaction is a natural generalization of the spin one half one-photon potential (see, 
e.g., (l3|). The only difference is that like in the Lagrangian in Eq. ([2]) the coefficients 
in Eq. (jl]) before the Darwin terms depend on the magnitude of particles' spins, and new 
terms with electric quadrupole moments arise. The interaction in Eq. (jlj) is calculated in 
the absence of the external magnetic field that is present in the f^-factor problem. This 
drawback is easily repaired by the minimal substitution pj — )■ pj — CjAj, Ai = B x rj/2. 

Combining the nonrecoil Lagrangian in Eq. ([2]) and the one-photon potential (after min- 
imal substitution) in Eq. (^ we obtain a total effective two-particle NRQED Hamiltonian 
for electromagnetically interacting particles with arbitrary spins (we preserve only the terms 
relevant for calculation of the g-iactor contributions) 



where 
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and H2 is obtained from H^ by the substitution 1 — )■ 2. 

The nonrelativistic effective two-particle Hamiltonian describes all (nonrecoil and recoil) 
leading relativistic corrections to bound state (/-factors of each of the constituents. To 
calculate these corrections we need to separate effects of the bound system motion as a 
whole from the internal effects. This task is not quite trivial because the center of mass 
variables do not separate in the presence of external field. For the current case of a small 
magnetic field a solution was suggested in p, [jU|. The main idea is to insist that the center 
mass of a loosely bound system moves in an external field exactly in the same way as a 
respective elementary particle with the same mass and charge. To satisfy this transparent 
physical requirement transition to the standard center of mass coordinates r = Vi — r2, 
R = /xiT"! + /i2^2, f^i = 'fni/ijni +^712) should be accompanied by the unitary transformation 
of the Hamiltonian H — )■ U^^HU, where U = e*('^i'^2-e2/ii)A(fl)r y^f^^gj^ i^j^jg transformation 
we extract from the transformed Hamiltonian terms that describe the spin interaction with 
an external field. The final Hamiltonian for the first particle is 
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The Hamiltonian for the second constituent has the same form. The leading binding correc- 
tion to the g-iactoT is completely described by this Hamiltonian. We calculate its matrix ele- 
ment with the help of the first order perturbation theory between the Schrodinger-Coulomb 
wave functions that are eigenfunctions of the unperturbed internal Hamiltonian. After sim- 
ple calculations we obtain the bound state (^-factors with account of the leading relativistic 
corrections of order {Za)"^ for s-states with the principal quantum number n 
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These results resolve the discrepancy mentioned in the Introduction in favor of the results 
in (6| (see also 15|)- The remarkable property of the expressions in Eq. (^ and Eq. (^ is 
that they are universal for particles of any spin; they depend only on the (/-factors of free 
charged particle, not on the magnitude of their spin. Technically this happened because all 
terms in the effective two-particle NRQED Hamiltonian in Eq. ([5]) relevant for calculation of 
the leading relativistic corrections do not contain spin-dependent coefficients Aj, Sf. On the 
other hand analysis of dimensions and spin structure of all terms in the NRQED Lagrangian 
in Eq. 02]) leads to the conclusion that terms with derivatives of electric fields do not generate 
contributions to the leading relativistic corrections to the bound state gf-factors. Omission 
of the field derivatives is the basic assumption for validity of the Bargmann-Michel-Telegdi 



(BMT) equation [13|, |16[. Hence, the approximation based on the BMT equation [61] is 
sufficient for calculation of the the leading nonrecoil relativistic corrections to bound state 
g-iactoTS. Then the leading relativistic corrections are universal because the BMT equation 
is universal for all spins. For purely electromagnetically interacting particles the free g- 
factors in the BMT equation and in Eqs. (iQl fTOj) are just sums of the QED perturbation 
series. The BMT equation alone is insufficient for calculation of the leading relativistic 
recoil corrections and should be amended by the one-photon exchange potential in Eq. ([7]). 
This is a spin-orbit interaction, it does not depend on magnitude of spins, and it produces 
universal corrections to bound state (/-factors. 
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